Solvation Effects in Near-Critical Binary Mixtures 
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A Ginzburg-Landau theory is presented to investigate solvation effects in near-critical polar fluid 
binary mixtures. Concentration-dependence of the dielectric constant gives rise to a shell region 
around a charged particle within which solvation occurs preferentially. As the critical point is 
approached, the concentration has a long-range Ornstein-Zernike tail representing strong critical 
electrostriction. If salt is added, strong coupling arises among the critical fluctuations and the ions. 
The structure factors of the critical fluctuations and the charge density are calculated and the phase 
transition behavior is discussed. 



I. INTRODUCTION 

Solvation effects are of great importance in under- 
standing the degree of solubility of ions in various polar 
solvents [1,2]. A large number of papers have been de- 
voted on this problem [3-6]. The theoretical approaches 
range from solving phcnomenological continuum models 
to performing computer simulations on microscopic mod- 
els. Originally, Born used a simple continuum theory of 
electrostatics to derive the solvation (polarization) free 
energy of a single ion with charge Ze [1,7], 



AG B = -{Z 2 e 2 /2R B ){\ - 1/e), 



(1.1) 



where Rb represents the ionic radius [1,8]. In his theory 
the solvent dielectric constant e is assumed to be homoge- 
neous and the contribution without polarization (e = 1) 
or in vacuum is subtracted. In the line of continuum 
models, considerable efforts have been made to take into 
account the dielectric saturation near the ion core (which 
are particularly crucial for multivalent ions) [9,10]. The 
Born theory also neglects possible inhomogeneity of s in 
the vicinity of the ion due to a change in the density for 
one-component fluids or in the concentration for binary 
fluid mixtures (electrostriction) [11-14]. 

We point out that the previous theories have not yet 
treated the physics on mesoscopic scales such as solva- 
tion effects on phase transitions or collective phenomena 
at non-dilute salt concentrations. In this paper we will 
show that such effects can be conveniently studied within 
a scheme of Ginzburg-Landau theory. First, we will ap- 
ply this approach to solvation around a charged particle 
in polar binary mixtures near the consolute critical point. 
Here important is the gradient free energy well-known in 
the field of critical phenomena [15], which will reasonably 
accounts for the long-range electrostriction around ions. 
In the previous theories [11-13] in the physical chemistry, 
the gradient free energy has been neglected. On the basis 
of these results we will then discuss the effects of low- 
density salt in polar binary mixtures in a Debye-Huckel 
type approximation. This second part is an extension of 
a previous paper by one of the present authors [16]. We 
also note that similar effects have recently been discussed 
when the solvent is a liquid crystal [17], where the direc- 



tor field can be much deformed around a charged particle 
over a long distance. 



II. SOLVATION AROUND A SINGLE CHARGED 
PARTICLE 

A. Electrostatics 

We place a charged particle in a near-critical binary 
mixture in a one-phase state in equilibrium, in which the 
dielectric constant e = e(<j>) depends strongly on the con- 
centration. This is the case when the two components 
have very different dielectric constants, ea and Eb, with 
ea considerably larger than eb- There is no established 
theory of the dielectric constant of such mixtures. How- 
ever, empirically, e can be expressed roughly as a linear 
function of the concentration in many relevant mixtures 
[2]. In the present theory, for simplicity, we assume the 
linear dependence, 



£ = £q + eu 



(2.1) 



where Eq and E\ are constants (with eb = £o an d 
Ea = £q + £i) and <f> is the concentration or the vol- 
ume fraction of the preferred component A. Debye and 
Kleboth [18] performed a light scattering experiment on 
a mixture of nitrobenzene (NB)+ 2,2,4-trimethylpentane 
in electric field, for which e(0) increased from 2.1 to 34.2 
with increasing NB with d 2 £/d4> 2 = 28.7 near the criti- 
cal point. The measured curve of e(4>) vs <j> may be fairly 
fitted to the linear form (2.1) (and fitting is better in 
the region <f> >4> c ). In their experiment the critical tem- 
perature was 29.16 C and the critical concentration was 
51.7 wt% of NB (corresponding to a volume fraction of 
0.381 of NB). In this paper the dielectric constant is as- 
sumed to be independent of electric field, although this 
is a questionable assumption in the vicinity of small ions 
[9,10]. 

The electric field E = — is induced by the electric 
charges and the electric potential $ satisfies 



-47rp(r) 



where p(r) is the charge density, 
space-dependence of £ is crucial. 



(2.2) 

In our problem the 
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B. Ginzburg-Landau Free Energy 

We construct a simple theory for a nearly incompress- 
ible mixture near the consolute critical point at a given 
pressure. Then ip = <p — <p c is the order parameter of 
liquid-liquid phase transition, where </> c is the critical con- 
centration. The total Ginzburg-Landau free energy of the 
system is given by 



/ 



dr 



f 



C 



— E 2 

8tt 



(2.3) 



The molecules of the two components have a common vol- 
ume vq = a 3 , although they can have different volumes 
in general. The free energy density / = f(<j>) is given by 
the Hildcbrand [19] (or Bragg- Williams [15]) expression, 



/ 



k B T 



.ln^+(l-$ln(l-$+x#l-$ 



> (2-4) 



where \ is the interaction parameter dependent on the 
temperature T. The critical-point values of <p and \ are 



1/2, Xc = 2, 



(2.5) 



respectively, in the mean field theory. The order param- 
eter is defined by 



ip = 0-1/2. 



(2.6) 



The coefficient C of the gradient term is of order k B T/a 
and we will use a simple form [20,21], 



C = k B T X /a. 



(2.7) 



The last term in F is the electrostatic contribution. The 
free energy functional (2.3) supplemented with (2.2) can 
be used generally in the presence of ions in polar near- 
critical fluids, 



C. Preferential solvation 

We consider a single charged particle at the origin of 
the reference frame. The charged particle has a radius 
R. It is convenient to assume that the charge density is 
homogeneous within the sphere as 



p = Ze/(ATrR 3 /3) (r < R) 



(2.8) 



and vanishes outside the sphere r > R. The total charge 
is given by Ze. In our theory R is a phcnomcnological 
parameter. Then the volume fraction <j> and the electric 
potential $ can be defined even within the sphere and, for 
the single charge case, all the quantities may be assumed 
to depend on space as 



(j> = <f>(r), $ = $(r), E = -$'(r)r _1 r, 



(2.9) 



where r = |r| and $' = d&/dr. The concentration <fi(r) 
decrease from 0(0) to 0oo with increasing r (see Fig.l to 
be presented below). In this case (2.2) is solved to give 

-e{4>)&{r) = ^ f dr ir 2 p(n) = §0(r), (2.10) 
~ Jo r 

where 9{r) = 1 for r > R and 

6(r) = {r/Rf (r < R) (2.11) 

under the assumption (2.8). The electrostatic free energy 
is now simply of the form, 



/ dr—eE 2 
J 8tt 



1 



1 



Z 2 e 2 / dr—9(r) 



(2.12) 



If e is replaced by the dielectric constant far from the 
ion Eoo, the above quantity is integrated to become a 
Born contribution (see (1.1)) [7], 



Fboo — ■^ 2 e 2 /2e 00 i?B- 



(2.13) 



In terms of the critical-point dielectric constant e c = 
£o + £i/2 and the order parameter t/'oo = <Poo — 1/2 far 
from the ion, we have 



£oo — £o + £l4>oo — £c + £lipc 
The Born radius is given by 

R B = 5i?/6. 



(2.14) 



(2.15) 



In the integral (2.12) the contribution from the ion inte- 
rior r < R is 1/5 of that from the ion exterior, resulting 
in the factor 1/(1 + 1/5) = 5/6 in (2.15). However, if we 
replace e by its value at the ion e(0) = e + £i</>(0), we 
obtain another estimate, 



F B0 = Z 2 e 2 /2e(0)R B 

= F Boo e c /[e + £i0(O)]. 



(2.16) 



As will be calculated in the appendix, the equilibrium 
free energy of a single ion accounting for the inhomoge- 
neous <j) and e in the mean field theory reads 



F sol = 4tt / drr 2 [f + C{<P') 2 /2 + sE 2 /8w} 



I 



dr 



Cr 2 ((b') 2 + 



2 e 2 d 



8-Kr 4 dr 



where 



/W = /(0)-/to»)-M0-&o), 



(2.17) 



(2.18) 



with fj, — /'(</>oo) being the chemical potential difference. 
As shown in the second line of (2.17), F so i consists of the 
gradient contribution (oc C) and the electrostatic contri- 
bution (cx e 2 ). The latter becomes the Born contribution 
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(2.13) or (2.16) if s is replaced by Soo or e(0). In Fig. 7 
below we shall see the relation Fbo < F so i < Fboo • 

With the aid of (2.11) the minimization condition 
SF/S(f> = n =const. in equilibrium is rewritten as 



In 



1 + 2^ 



- 2\ip - xa 2 V 2 ip - v = A 



4/)2 



a 4 9 



\-2ip_ 

where e is the normalized dielectric constant, 
i = e/e c = 1 + (ei/e c )ip. 



(2.19) 



(2.20) 



The dimensionless chemical potential difference v 
von/k^T is written in terms of V'oo = tp{oo) as 



v = ln[(l + 2V>oo)/(l - 2^)] - 2x^oo. 



(2.21) 



The coefficient A on the right-hand side of (2.19) repre- 
sents the strength of solvation defined by 



A = e 1 Z 2 e 2 /8irs 2 c k B Ta 
= eiZ 2 £B/STre c a. 

In the second line the Bjerrum length 

£ B = e 2 /e c k B T 



(2.22) 
(2.23) 



at the critical concentration is introduced. For exam- 
ple, if we set £b ~ 14A, a ~ 2A, ei/eo ~ 2, we obtain 
A ~ Z 2 /5. In the vicinity of the ion the right-hand side 
of (2.19) takes a maximum of order Aa 4 /R 4 (because 
i <~ 1). If this maximum is much larger than 1, there 
arises a solvation shell around the ion within which pref- 
erential solvation is strong and <p = 1 . Thus the condition 
of strong solvation may be expressed as 



A > {R/af. 



(2.24) 



The solvation is stronger for smaller R/a, which is a well- 
known result in the literature [1]. 

It is worth noting that Padova presented an equation 
similar to (2.19) for the concentration in the vicinity of 
an ion using thermodynamic arguments [11], where C = 
asnd the gradient term (oc V 2 V>) did appear. Even far 
from the critical point, as can be seen in Fig. 8 below, the 
result with C — is considerably different from that with 
the gradient term in the strong solvation case. 



D. Relations near the critical point 

If the system is near the critical point far from the 
charged particle, we have \ — 2 and (p^ = 1/2. Then 
the inverse correlation length k = far from the ion is 
determined by 



aV = 4/(1 - 4</4) X - 2 



S2-y + 8$o. 



(2.25) 



As T — > T c in the UCST case we may set 
2-x = #i(T/T c -l). 



(2.26) 



The correlation length in the mean-field expression be- 
comes £ = i 0+ (T/T c - l)- 1 ' 2 with 



£o+ = aD 



-1/2 



(2.27) 



in one-phase states at the critical composition [18]. The 
fluid can remain near-critical only in the region where 
Aa 4 /r 4 < 1 or 



r > A l/i a. 



(2.28) 



The shell radius should thus grow as A 1 / 4 a for A S> 1. 
In the near-critical region, (2.19) may be approximated 
in the Landau expansion form [16], 



a 2 [n 2 - V 2 ]<5V + o( 3 ^oo^ 2 + H 3 ) = A 



a 4 9 2 
2r 4 ' 



(2.29) 



where Sip — xjj — ip^ . This equation does not hold within 
the solvation shell if ip is not small there for strong sol- 
vation, but holds generally far from the ion. If we retain 
the linear terms only in the above equation, we obtain 
(k 2 — V 2 )<5?/> = far from the ion, so there should emerge 
the Ornstcin-Zernike tail, 



8ip(r) ^ — e 
r 



(2.30) 



in the concentration. The coefficient B is a microscopic 
length to be determined in the following. 

In the Landau theory of phase transition (usually pre- 
sented in the context of Ising spin systems) [15] , the right- 
hand side of (2.19) or (2.29) plays the role of a (position- 
dependent) magnetic field h. If h were homogeneous and 
small, we would have the linear response Sip ~ h/n 2 . 
In the present case, however, h = h(r) steeply changes 
in space (oc r~ 4 ) from a maximum of order yla 4 /i? 4 to 
zero and as a result the gradient term — V 2 ip becomes 
increasingly crucial with decreasing k. 



E. Weak solvation limit near the critical point 

It is instructive first to consider the simplest case of 
very small A near the critical point where < 2 — \ <C 1 
and Voo = 0. Furthermore, if ip(0) <C na holds at the 
ion, we may neglect the second term on the left-hand 
side of (2.29) in the whole space to obtain the linearized 
equation, 



(k 2 - V 2 )iP = A9(r) 2 a 2 /2r 4 



(2.31) 



where we have set i = 1. This equation is readily inte- 
grated to give 
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ip(r) = 



Aa 2 
4nr 



dr-i 



e- Klr - ril 6{n) 



D 



(2.32) 



where the coefficient D is determined such that ^(0) is 
finite. We can see the relations, 



V>(0) = 2kD = Aa 2 / dr 1 e- Kri 9(r 1 ) 2 /2rf 
Jo 

= 3Aa 2 /8R 2 (kR^I). (2.33) 

The second line holds for kR <C 1. As stated above 
(2.31), the linearization is allowed for 



F. Behavior at the critical point 

At the critical point (x = 2 and = 0.5), i[)(r) is of 
the form Ba/r at long distances in the linear theory. We 
here show that it is delicately modified in the nonlinear 
theory. We rewrite (2.29) for r 3> a as 



(2.40) 



A < kR 2 /a or eiZ 2 £ B /8Tre c < R 2 n. 



(2.34) 



We set G(s) — rtp(r)/a with s = ln(l + r/a) to obtain 



From (2.31) the three-dimensional Fourier transforma- 
tion of ip becomes Vfe = H k /(n 2 + k 2 ) where 



(i-O 2 



tf__d_ 

ds 2 ds 



G( S ) = -G( S Y--Ae-s. (2.41) 



H k = 2nAa 2 / drd(r) 2 sm(kr)/kr 3 . (2.35) For s > 1 we assume the algebraic decay G(s) 

Jo s~@. Then the above equation is approximated as 



Here H k tends to a well-defined long-wavelength limit 
H = limfc^o-fffc for kR < 1. Thus, for kR < 1, ip(r) 
approximately takes the Ornstein-Zernike form (2.30) in 
a wide space region with B = a 2 _ff /47r. Here, 



B = Aa 2 



rOC 

/ dr9(r) 2 /2r 2 
Jo 



Aa 2 /2R B , 



(2.36) 



with R B being given by (2.15). This form also directly 
follows from the integral form (2.32) (by taking the limit 
k — ► 0) . If we do not assume the special form C in 
(2.7), we have B = {k B T/C)Aa/R B (which will be used 
in (3.12) below). 

At very long distances r > A, however, the gradient 
term — V 2 ?/> in (2.31) becomes negligible to give 



ip{r) S Aa 2 /2n 2 r 4 . 



(2.37) 



This is the relation of linear response as stated at the 
end of subsection 3D. The crossover length A is longer 
than £ = k^ 1 and is determined by the balance of the 
Ornstein-Zernike form (2.30) and the long-distance limit 
(2.37) as 



(«A) d e 



3 -kA 



= kR. 



(2.38) 



For example, A = 10. 8£ for R = 0.03£. Thus kX increases 
as kR decreases. 

We may then calculate the excess accumulation ^ of 
the component A around the ion as 



* = J drtp(r) 

= 2nAa 2 /K 2 R B = AttB/k 2 , 



(2.39) 



which coincides with the space integral of the Ornstein- 
Zernike form (2.30) and grows strongly on approaching 
the critical point. The integral of the tail (2.38) in the 
region r > A is much smaller than the value in (2.39) by 
R/X and is canceled with a contribution from the correc- 
tion to the Ornstein-Zernike form in the range r < X. 



-dG(s)/ds = 8G(s) 3 /3, yielding the solution 
G(s) £* {&l 2 /A)s- l l 2 



(2.42) 



with ft = 1/2. Thus the asymptotic behavior at ex- 
tremely long distances s >• 1 should be 



i/)(r) = 3 1/2 a/[4ry/\n(l + r/a)}. 



(2.43) 



The above form will turn out to be consistent with nu- 
merical analysis in Fig. 6 to be presented below. 



G. Numerical Results 



We then show numerical solutions of (2.19) around an 
ion in equilibrium one-phase states for the case e = 10 
and £i = 70 in (2.1), assuming the model charge distribu- 
tion (2.8). The ion radius is given by R = 0.3a in terms 
of the solvent radius a. 

We show the concentration profile, <p(r) vs r/a, in Fig.l 
and the concentration 0(0) at the ion center in Fig. 2. No- 
tice that (j>(r) = (j>(0) within the interior of the ion r <R 
owing to the smoothing effect of the gradient term. We 
can see that the concentration <p(0) at the ion increases 
obeying the linear growth (2.31) for A <C 1 and saturates 
into 1 for A >0.3 and that the shell region expands with 
increasing A. The shell radius i? s hcii may be defined such 
that 4>{r) > 0.75 for r < R s hdh sa Y- Then i? s hcii ~ A a 
with the effective exponent a being about 1/4 for the 
profiles (d)~(g) in Fig.l. Below (2.28) we have discussed 
that a should tend to 1/4 for A > 1. 

In Fig. 3 we examine the behavior of ip(r) = <p(r) — 1/2 
with = 1/2 at A = 0.346. We can see that the com- 
bination rip(r)/a behaves as (B/a) exp(— nr) in a wide 
intermediate range a <r < A. The correlation length 
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£ = k 1 is given by 10a for \ = 1-99 and by 10 3 / 2 a for 
X = 1.999. For r >A, V( r ) can be fitted to the r~ 4 tail in 
(2.37). Fig.4 displays B/a vs A for x = 1-99 and 1.999 at 
0oo = 0.5. In the nonlinear regime A >0.1 the data may 

be fitted to B/a ~ A - 4 . In Fig. 5 we plot the crossover 
length A divided by a obtained from equating the two 
fitted functions (2.30) and (2.37). In the linear regime 
with small A, A is given by (2.37). 

The concentration profile just at the critical point is 
also of great interest. In Fig. 6 we plot the combination 
V'(r)rs 1 / 2 vs s = ln(l + r/a) in the region s < 12 or 



r/a < 



1 for A = 0.2,1, and 5. If (2.42) is valid, 



this combination should tend to 3 1 / 2 /4 = 0.43. Thus the 
three curves are consistent with (2.42). 

Fig. 7 displays the numerically calculated solvation free 
energy F so i in (2.17) divided by F Boo in (2.13). Also plot- 
ted is the ratio -Fbo/^Boo = e c /[eo + £10(0)] in (2.26). 
Thus we find that F so i is intermediate between Fbo and 
Fboo- Coincidence of the two curves of x = 1-99 and 
1.999 implies that F so i tends to a constant on approach- 
ing the critical point. 

In Fig. 8 we demonstrate that rip(r)/a vs r/a far above 
the critical point with x = 1-5 and (f)^ = 0.5. The pro- 
file can be fitted neither to the Ornstein-Zcrnike form 
nor to the solution of (2.19) without the gradient term. 
The curves from the local equilibrium theory without the 
gradient term [11] considerably differs from that in our 
theory even away from the critical point. 



III. NEAR-CRITICAL FLUIDS WITH IONS 



Experimentally, it has long been known that even a 
small fraction of ions (salt) with c <C 1 dramatically 
changes the liquid-liquid phase behavior, where c is the 
mass or mole fraction of ions. For small c, the UCST 
coexistence curve shifts upward as 



AT ion = A ion c + 0(c 2 ) 



with large positive coefficient Ai on , expanding the region 
of demixing. For example, Ai on /T c ~ 10 with T c ~ 300K 
when NaCl was added to cyclohcxane+mcthyl alcohol 
[22] and to triethylamine+H 2 [23]. The LCST co- 
existence curve of 2,6-lutidine+H20(D20), the shift is 
downward with |Ai on |/T c ~ 10 [24]. Similar large im- 
purity effects were observed when water was added to 
methanol+cyclohexane [25]. In some aqueous mixtures, 
even if they are miscible at all T at atmosphere pressure 
without salt, addition of a small amount of salt gives 
rise to reentrant phase separation behavior [26-29] . Such 
reentrant phase behavior is believed to arise from hydro- 
gen bonding. 



A. Ginzburg-Landau free energy 

In the previous section we have examined the concen- 
tration profile around a single fixed charged particle. We 
then need to construct a theory on polar binary mix- 
tures in which a small amount of salt is doped [16]. Even 
if the ion concentration is very small, the interactions 
among the ions should be much influenced by the prefer- 
ential solvation near the critical point. We here present 
a Ginzburg-Landau free energy in which the fluctuations 
on the spatial scale of the solvation shell radius i? s h e ii 
have been coarse-grained. Then the order parameter 
ip(r) = 4>{r) — 4> c and the ion densities ni(r) and n 2 (r) 
for the two species do not change appreciably on the scale 

Of -Rshcll- 

The two ion species, 1 and 2, have charges, Q\ = Ze 
and Q2 = — e. The average densities are written as 



(n 2 ) = Zn, 



(3.2) 



where (• • •) denotes taking the space average. Notice that 
the charge neutrality condition yields Z(m) = (n 2 ). The 
total charge number density is written as 



n to t = (Z + l)n. 
We propose the free energy F in the form [16], 



(3.3) 



dr 



K=l,2 



(3.4) 



The first line is of the same form as the free energy in 
(2.3), but f{4>) can be expressed in the Landau expan- 
sion form / = rijj 2 /2 + uoip 4 /4 because the non-critical 
behavior in the solvation shell is not treated here. In 
addition we do not assume the special form (2.7) for the 
coefficient C. The electric field E = — here arises 
from the Coulomb interaction among the ions so that 



(3.1) where 



V • E = -V 2 $ = 4 7 r£- 1 / o(r), 



p(r) = e[Zni(r) - n 2 (r)] 



(3.5) 



(3.6) 



is the charge density. Therefore, the electrostatic part of 
the free energy is rewritten as 

Because the small-scale electrostriction on the scale of 
i? s hdi has been coarse-grained, we neglect the small-scale 
inhomogeneity of the dielectric constant and set e = e c 
with e c being the dielectric constant at the critical point. 

The free energy with the bilinear coupling (oc nxip) 
was already examined in Refs. [30,31], where the cou- 
pling constants wk were phcnomenological parameters, 
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however. We claim that the small-scale solvation or elec- 
trostriction is the origin of the bilinear coupling. In fact, 
the minimization of F at fixed inhomogeneous ion densi- 
ties yields 

/i = /' - CV 2 ip + k B T{w x ni + w 2 n 2 ) = const., (3.8) 

where /' = df /dip. For small inhomogeneous deviations 
Sip — ip — (V') aim Stik = n-K — (tik) (K = 1,2) we 
linearize the above equation to obtain 

C(k 2 - V 2 )5ijj + fc B T(wi<5ni + w 2 5n 2 ) = 0, (3.9) 

where Ck 2 = d 2 f /dip 2 with k being the inverse correla- 
tion length. Notice that, if the ion densities are localized 
in regions shorter than £ = k _1 , Sip has Ornstein-Zernike 
tails in the form of (2.30) around such localized regions. 
Thus wk is written as 



w K = -(4nC/k B T)B K (K = 1, 2), 



(3.10) 



where Bk is the coefficient of the Ornstein-Zernike tail 
around an ion of the species K. If the two components 
of the fluid mixture have the same molecular size a as 
assumed in the previous section, Bk is determined by 
the parameters Ak in (2.22), where 



Ax/Z 2 = A 2 = ei^ B /87T£ c a. 



(3.11) 



In particular, in the linear solvation regime with small 
A K , (2.36) yields 



wk = — 47rAKa/ Rbk (linear), 



(3.12) 



where Rbk is the Born radius of the ion species K. 
See Fig. 4 for the nonlinear behavior of Bk- That is, 
Bk /a <~ A K 4 for large A k and % can well exceed 1 . 

In Ref. [16] wk was related to the Born solvation free 
energy AG BK = -(Q 2 K /2R BK ){1 - l/e) in (1.1) which 
should be defined for the two species. That is, treating e 
as a function of ip, we obtain 

Q 

k B Tw K = ^AG Bif - - £l Q 2 K/2e 2 c R BK , (3.13) 

where the differentiation is taken at <f> = (p c . Obviously, 
(3.12) and (3.13) are equivalent from (2.22). This simple 
derivation cannot be used in the strong solvation regime, 
however. 



B. Fluctuations in one-phase states 

We consider small plane-wave fluctuations in a one- 
phase state. The fluctuation contributions to F in the 
bilinear order are written as 



SF 



\ { r + Cf)\ip q ? + ^- 2 \ Pq ? 



(3.14) 



where J" • • • = (2n) 3 J dq ■ ■ ■ denotes the integration 
over the wave vector q and 



r = d 2 f/diP 2 . 



(3.15) 



The ipq, nnq, and pq are the Fourier transformations of 
Sip, SriK, and p, respectively. We may set 



r = a Q (T-T c 



(3.16) 



at the critical concentration, where the coefficient ao 
is of order k B /a 3 and is equal to Ak B Di/a 3 with Di 
being defined by (2.26). In the mean-field expression 
£ = £o+(T/T c - I)- 1 / 2 of the correlation length, we have 



a = C/^ 2 + T c . 



(3.17) 



We may then calculate the structure factor for the 
order parameter S(q) = {\tpq\ 2 ) in the mean field the- 
ory. Since the equilibrium distribution is given by const. 
cxp(— SF/K B T), the inverse of S(q) is written as [30,16] 



k B T/S(q)=a (T-T c -AT ion ) 

+Cq 2 [l- 1 2 J(l + \lq 2 ) 



(3.18) 



where ATj on is the critical temperature shift due to salt 
given by 

AT ion = {k B T c /a )( Wl + Zw 2 ) 2 n tot /(l + Z) 2 , (3.19) 
and Ad is the Debye screening length defined by 

Ad 1 = ^Zn tot e 2 /e c k B Tfl 2 . 
We introduce a dimensionless parameter, 



(3.20) 



7p = K - w 2 \(k B T/AnCe B ) 1/2 /(l + Z) 
= \B X - B 2 \(4nC/£ B k B T) 1 / 2 /(l + Z). 



(3.21) 



This number is independent of the ion density and rep- 
resents the strength of asymmetry in the ion-induced po- 
larization between the two components. In the linear 
solvation regime with small Ak, we use (3.12) to obtain 



7 P = 



2e c (l + Z)VI^C 



Z 2 
Rbi 



1 

RB2 



(linear). (3.22) 



The structure factor of the charge density fluctuations 
Spp(q) = {\pq\ 2 )/e 2 is written as 



S pp(i) = Zn t 



+(wi - w 2 y 



A D g 2 
A + \ 2 D q 2 



Zn t . 



X 2 n 2 



S(q), (3.23) 



6 



where the first term is the Debye-Hiickel structure factor 
and the second term arises from the coupling to the order 
parameter. 

From the structure factor S(q) in (3.18) we may draw 
the following conclusions. 

(i) If 7 P < 1, S(q) is maximum at q = and the criti- 
cal temperature shift due to ions is given by (3.19). As 
a rough estimate for the monovalent case Z = 1, we 
set a^n ~ k^T^_ n ~ c, where c is the mass or mole 
fraction of the ions. Then ATi on <~ (wi + W2) 2 c. If 
\wi + w 2 \ ~ 3, this result is consistent with the experi- 
ments [22,23]. In future experiments, let k^T '/CS(q) vs 
q 2 be plotted; then, the slope is 1 — 7 p for qX^ <C 1 and 
is 1 for q\o 3> 1- This changeover should be detectable 
unless 7 P <C 1. 

(ii) The case 7 P = 1 corresponds to a Lifshitz point 
[32,33], where l/S{q) - 1/5(0) cx q 4 /(l + q 2 \ 2 D ). 

(iii) If 7 P > 1, the structure factor attains a maximum at 
an intermediate wave number q m given by 

g m = (7p-i) 1/2 AD, (3.24) 

so q m cx ritot 1 ^ 2 - The maximum of the structure factor 
S(q m ) is written as 

S(q m ) = k B T/[a (T - T c - AT( on )}, (3.25) 

where 

AT/ on = AT ion + ( 7p - 1) 2 T C ^ + /A£, (3.26) 

where use has been made of (3.17). A charge-density- 
wave phase should be realized for T — T c < ATj' on . It is 
remarkable that this mesoscopic phase appears however 
small ntot is as long as 7 P > 1 and ij m L » 1 with L 
being the system length. Here relevant is the coupling of 
the order parameter and the charge density in the form 
cx tpp in the free energy density, which generally exists 
in ionic systems. This possibility of a mesoscopic phase 
was first predicted for electrolytes by Nabutovskii et al 
[30], but has not yet been confirmed in experiments. In 
polyclcctrolytes, on the other hand, electric charges are 
attached to polymers and the structure factor of the poly- 
mer is known to take a form similar to (3.18), leading to 
a mesoscopic phase at low temperatures, in the Debye- 
Hiickel approximation [32]. 

IV. SUMMARY AND CONCLUDING REMARKS 

We summarize our results, 
(i) In the first part, we have examined the concentra- 
tion profile around a charged particle in a near-critical 
polar binary mixture, as can be seen in Figs.l and 2. 
We have started with the simple Ginzburg-Landau free 
energy (2.3) supplemented with the relations (2.1) and 
(2.2) of electrostatics. Preferential solvation can occur in 



the presence of a concentration-dependent dielectric con- 
stant and becomes strong when the parameter A defined 
by (2.22) satisfies the condition (2.24). As the critical 
point is approached, the concentration has a long-range 
Ornstcin-Zernike tail (2.30) in an intermediate range 
a < r < A as demonstrated in Fig. 3. We plot the co- 
efficient B vs A in Fig. 4 and the crossover length A vs A 
in Fig. 5. At the critical point ip(r) is long-ranged in the 
form of (2.43) as demonstrated in Fig. 6. The solvation 
free energy F so \ behaves as a function of A as in Fig. 7. 
(ii) In the second part of this work we have presented the 
coarse-grained free energy as in (3.4) where the order pa- 
rameter and the ion densities are strongly coupled in the 
bilinear form. We have related the coupling constants 
wi and w 2 in terms of B\ and £> 2 , the coefficients of the 
Ornstcin-Zernike tails for the ion, 1 and 2, respectively. 
Furthermore, if this coupling is significantly different be- 
tween the two species of the ions, a charge-density-wave 
phase can be realized. 

We then give some remarks. 

(i) We have neglected the dielectric saturation near the 
ion [9,10]. If this effect is taken into account, the growth 
of the coefficient B with increasing A would be weaker 
than in Fig. 4. 

(ii) In future experiments, the structure factor (3.18) 
should be observed. There could be the case 7 P > 1, 
where the mesoscopic phase is formed at low temper- 
atures, for pairs of small cations like Li + or Al 3+ and 
relatively large anions [1] . 

(iii) Dynamical problems remain largely unsolved. For 
example, we are interested in the response of ions against 
ac electric field in near-critical polar binary mixtures. In 
the strong solvation condition the ion mobility should 
strongly depend on lot^ where lu is the frequency of 
the field and is the order parameter life time (— 
Qnr]^ 3 /kBT with r\ being the shear viscosity [15]). Fur- 
thermore, even if uit^ <C 1, the large-scale solvation cloud 
on the scale of £ should become nonlincarly dependent on 
the applied field. 

(iv) Solvation or polarization effects due to inhomogc- 
neous dielectric constant (or dielectric tensor for liq- 
uid crystals) can strongly influence the phase transition 
behavior in complex fluids including polyelectrolytcs, 
charged gels, and liquid crystals containing ions or 
charged colloids [17]. The importance of this effect in 
polyelectrolytes has recently been pointed out by Kra- 
marenko et al. [34] 

(v) A similar Ginzburg-Landau theory can easily be de- 
veloped on solvation effects in polar one-component flu- 
ids. For example, we may start with the free energy 
density in the van der Waals theory [15] in place of the 
Hildebrand free energy density (2.4). As regards statics 
salient features in one-component fluids are nearly the 
same as those in polar binary mixtures, but dynamics of 
critical electrostriction should be very different between 
the two cases [15]. 
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Appendix 

In equilibrium <j>(r) around a charged particle is deter- 
mined by 

/'-^-^W,0, (A.l) 

where / is defined in (2.16). Multiplying <f> = d<fr/dr and 
integrating over r in the range [r, oo] , we obtain 

jf^W. <A.2) 

where 

„. , 2C. „., Z 2 e 2 d [6 2 \ 

H(r) = — (4>) + 7- -r • ( A -3) 

W r KV ' 8ire dr \r 4 J V ; 

Integration of (A. 2) over the whole space yields (2.17) 
with the aid of f °° drr 2 dnff (n) = J °° drr 3 H(r)/3. 
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r/a 

FIG. 3. Numerical rtp(r)/a vs r/a on a semi- logarithmic scale at A = 0.346 for \ = 1.99 (lower solid line) and 1.999 
(upper solid line) with 0oo = 1/2. The numerical tp(r) = 4>(r) — 1/2 can be excellently fitted to the Ornstein-Zernike form 
(corresponding to dotted lines) in the intermediate range a <r <X, while it can be fitted to the r -4 tail (2.37) (corresponding 
to broken lines) for r >\. The curve of A vs A will be given in Fig.5. 



1 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — r 




-2-10 1 
log 10 A 

FIG. 4. Coefficient £> of the Ornstein-Zernike tail in the concentration profile as a function of A at (poo — 0.5 for \ — 1-99 
and 1.999. It exhibits only weak dependence on 2 - x « T - T c . 
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-3 -2-10 1 



log 10 A 

FIG. 5. Crossover length A from the Ornstein-Zernike tail (2.30) to the r~ 4 tail (2.37) at 0oo = 0.5 for \ ~ 1-99 where 

i = ioo. 
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